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ABSTRACT 


The  windowed  discrete  Fourier  transform  (DFT)  is  a  widely  used  tool  in  count¬ 
less  signal  processing  applications.  Surprisingly,  there  lias  been  little1  analysis  of  its 
accuracy  for  the  processing  of  random  bandliniited  signals.  This  report  derives  an 
error  formula  at  leading  asymptotic  order  in  the  number  of  data  samples.  The  for¬ 
mula  applies  to  essentially  any  window  function,  with  explicit  results  tabulated  for 
some  of  the  most  common  eases.  The  asymptotic  error  is  shown  to  agree  well  with 
Monte  Carlo  results,  even  for  very  small  numbers  of  samples. 
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1.  INTRODUCTION 


Digital  signal  processing  is  largely  founded  on  the  Nvquist  theorem .  which  tells  us  that  a 
bandlimited  signal  can  he  recovered  from  its  unaliased,  discrete  samples.  The  Nvquist  theorem 
makes  it  possible  for  a  calculation  based  on  discretely  sampled  data  to  exactly  match  one  based  on 
continuous  quantities.  A  good  example  is  the  Fourier  transform  of  a  bandlimited  signal  s(t): 

r  =  At  jr  *(»*A t)e-'2"iltnAl.  (1) 

7l=  —  OC' 

The  equality  is  exact  when  A /  <  A/nvq,  where  A/nVq  is  the  Nvquist  time  step. 

Unfortunately,  Equation  (1)  does  not  apply  to  practical  computations  since  it  uses  ail  infinite.' 
number  of  samples.  Finite-duration  signals  cannot  be  exactly  bandlimited,  so  there  must  be  some 
amount  of  aliasing.  A  common  approach  is  to  use  a  smooth  windowing  function  w(t)  to  select  a 
finite  time  interval  T  from  the  bandlimited  signal: 


OC 

%  A /  ^  s(nAt)w 

n—  —  o c 


— 'InifinAt 


(2) 


The  window  function  tc(t)  is  assumed  to  be  real,  non- negative,  and  strictly  zero  when  |t|  > 
1/2.  If  one  is  free  to  use  long  time  intervals  T  and  a  smooth  window'  function,  the  error  in  Equation 
(2)  can  easily  he  made  very  small.  On  the  other  hand,  in  certain  applications  it  is  likely  that  T  will 
he  modest  or  small,  in  which  case  the  error  may  be  significant.  One  example  of  such  an  application 
is  joint  time-frequency  analysis:  one  can  be  forced  to  used  short  time  windows  when  the  signal  has 
rapidly  changing  spectral  content  [1]. 

Another  example,  which  was  the  original  motivation  for  this  work,  is  the  development  of  an 
efficient  backprojection  image  formation  algorithm  for  synthetic  aperture  radar  [2.4].  The  algorithm 
achieves  efficiency  by  decomposing  the  synthetic  aperture  integration  into  a  nested  hierarchy  of 
synthetic  apertures.  Maximum  efficiency  is  achieved  when  the  bottom  layer  of  the  hierarchy  consist  s 
of  integration  over  synthetic  apertures  of  only  a  few  Nvquist  samples.  Thus,  algorithmic  efficiency 
forces  one  into  the  case  where  the  error  in  Equation  (2)  is  likely  to  be  significant. 

The  error  in  Equation  (2)  may  also  be  of  concern  in  applications  where  the  required  accuracy 
is  very'  stringent.  A  common  example  is  adaptive  array  processing  where  one  wishes  to  cancel  very 
strong  interference  sources  in  order  to  detect  a  weak  signal  of  interest.  The  error  incurred  in  the 
finite  sum  of  Equation  (2)  will  act  like  a  noise  source  which  scales  with  the  strength  of  the  strongest 
signals.  Successful  detection  will  require  that  the  error  be  smaller  than  the  ratio  of  the  weak  signal 
to  the  strong  interference. 

In  [4],  the  authors  consider  a  version  of  Equation  (2)  with  a  square'  window  function.  An 
exact  error  expression  is  computed  for  deterministic  signals  with  piecewise'  polynomial  structure. 
For  more  general  classes  of  deterministic  signals,  the  same  error  expression  is  show  n  to  hold  in  tin1 
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limit  whore  At  — > -  0  and  //  — *  oo  while  [iAt  is  held  fixed.  While  the  results  are  mathematically 
correct,  they  arc'  not.  useful  for  the  typical  signal  processing  scenarios  because  they  do  not  capture 
the  effects  of  a  finite  number  of  data  samples.  For  example,  the  error  is  exactly  zero  for  the  zero 
frequency  component . 

In  this  report,  the  signal  s(t)  is  treated  as  a  stochastic  quantity  and  arbitrary  window  functions 
w(t)  are  allowed.  The  mean  squared  error  in  Equation  (2)  is  computed  at  leading  asymptotic  order 
in  the  number  of  samples  N  —  T  /At.  Explicit  results  are  given  for  many  common  window  functions. 
The  error  estimates  are  shown  to  agree'  well  with  Monte  Carlo  results,  even  when  the  number  of 
samples  is  relatively  small. 
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2.  DEFINITION  OF  THE  PROBLEM 


Let  s(t)  be  a  random  signal  of  bandwidth  B,  with  Fourier  transform  S(v): 

»{t)  =  f°°  S(v)c1%u'Ldv  (3) 

•/  -OO 

S{u)  =  0  for  |t/|  >  j  (-1) 

The  signal  s (t)  is  assumed  to  be  stationary  with  zero  mean  and  spectral  density  f)(v): 


E[.s(t)]  =  0 

[B/ 2  t ? 

E  [.s(/).s(/')*]  =  /  p{y)cl7Tli'^  1  ^dvp(v)  —  0  for  \u\  >  — 
J-B/2  ~ 

The  windowed  Fourier  transform  is  given  by  the*  integral 


0) 


(B) 


where  the  unit-interval  window  function  w(r)  is  assumed  to  be  real,  non-negative,  and  strictly  zero 
when  |t|  >  The  frequency  //  is  assumed  to  be  restricted  to  the  interval  |//|  <  ly.  The  discrete 
approximation  to  the  above  integral  is  given  by 


OO 

In  =  A f  ^  s(tn)w 

71  —  —  DC 


e-27T  ifltn 


(7) 


where  tn  =  (rr  +  a)  At  and  the  timo  step  At  meets  the  Nyqnist  condition:  At  <  A/„y,i  =  7j.  The 
parameter  o  controls  the  offset,  of  the  window  edges  with  respect  to  the  sample  locations  of  the 
discrete  data.  The  signal  is  not  assumed  to  be  synchronized  with  the  sampling  in  any  special  way. 
so  a  is  treated  as  a  random  variable  with  uniform  distribution  in  the  interval  [—  t-},  ^].  It  is  not 
necessary  to  consider  other  values  of  n  due  to  the  stationarity  of  s(t). 

Our  goal  is  to  estimate  the  normalized  mean  square  error,  defined  as 


o 

e 


E 


\I-In\ 2 


(«) 
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3.  EVALUATION  OF  EXPECTATION  VALUES 


Using  the  Poisson  summation  formula  and  other  standard  Fourier  identities,  I  and  //>  can  he 
rewritten  in  the  forms: 


=tF 


S{u)W  (T{fj  -  v))dv 


(!)) 


and 


OC  n 

h)  =  T  E  /  5(i/)^(t(/i-i/+^ 


ldi', 


(10) 


where  W(v)  is  the  Fourier  transform  of  the  window  function  </’(r).  losing  the  above  equations  and 
Equation  (5),  the  desired  expectation  values  may  he  computed: 


OC  n  'TXL, 

E  [|J  -  /op]  =  T2  Y,  J 


du 


1 1  —  CKj 


(11) 


and 


|/|-1  =  T2  ["  p{u)\W{T{tl-u))Fdu. 

L  J  J  -OC 


(12) 


!> 


4.  LEADING  ORDER  ASYMPTOTICS 


Up  to  this  point,  all  results  have  been  exact.  The  focus  will  now  shift  to  the  leading-order 
asymptotic  behavior  of  Expiation  (8)  when  the  window  size  T  becomes  largo  while  the  sampling 
spacing  At  and  frequency  //  are  held  constant. 

Consider  the  denominator  of  Equation  (8).  Making  the  change  of  variables  0  =  T(p  —  v)  in 
Equation  (12)  gives 


M3 


—  OO 
+  > 


T 


Tp(/'+)  I  \W(0)\2<W  +  Tp(,r)  r \W(O)\2d0 

J-oc  Jl) 


as  T  — *  oc.  The  notations  p(p^)  and  p(fj  )  are  to  be  understood  as  the  limit  of  p(i')-  v  — *  // 
approaching  from  above  and  below,  respectively.  Since  w(t)  is  strictly  real-valued,  the  Fourier 
transform  has  the  symmetry  |VT(0)|2  =  |Wr(— 6l)p.  The  leading-order  result  can  therefore  lx4  written 


u'(r)2c/r.  (13) 

•I  —oc 

For  the  numerator  of  Equal  ion  (8),  a  different  approach  is  necessary:  an  asymptotic  expansion 
will  be  derived  for  Vlr(/z),  then  plugged  in  to  Expiation  (11).  Note  that  w(r)  cannot  bo  analytic  for 
all  r  since  it  is  identically  zero  for  |r|  >  It  will  be  assumed  that  ir(r)  is  of  order  p  >  0,  meaning 
that: 

1.  ie^(r)  is  differentiable  everywhere  for  k  <  />; 

2.  w^(t)  is  differentiable  everywhere,  except  possibly  for  a  finite  set  of  point s  —  \  —  r \  <  r<>  < 
•  •  •  <  T/,  =  2* 

3.  vM\t)  must  lx*  discontinuous  at  at  least  one4  of  the  7>,  but  not  necessarily  all  of  them. 

It  would  be  difficult  to  imagine  any  useful  window  function  which  does  not  meet  the  above  condi¬ 
tions. 


E 


VY 


p(n+)  + 


Next,  write  W(v)  in  the  form 

W(v)  =  ^  I  w(r)c~27riTdr. 
k—  1 

Performing  integration  by  parts  repeatedly  and  collecting  terms  yields 

VV»  =  /n  *  y  w^(Tk)  +  - - T—  y  [  ’  (.  -2*.rr  tp+l)( 

v;  (2nus)i,+  l  "  (27 r»VF+1^yr, 


(14) 


where  A w^(r)  =  —  w^'\t  ).  The  Rieinann-Lebesgue  lemma  [5]  guarantees  that  the 

integrals  on  the  right-hand  side  vanish  as  u  —*  oc.  so  to  leading  order 


U» 


1 

(27t*V)p+* 


k=\ 


(15) 


Inserting  the  asymptotic  expression  above  into  Ec|uation  (11)  and  keeping  only  terms  of 
leading  order  as  T  — *  oo  yields: 


E 


\I-Id I2 


f-  r  pM  ,, 

(2 J-i-.  (/■  -  "  +  &P'+2  ’ 

rt/O 


(16) 


Note  that  in  expanding  |IE|2  via  Equation  (15).  all  terms  where  the  complex  exponentials  do  not 
cancel  have  been  dropped.  It  can  be  shown  via  the  Riemann-Lesbegue  lemma  that  the  dropped 
terms  are  of  subleading  order  after  integration  over  v. 
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RESULTS  FOR  FLAT  SPECTRAL  DENSITY 


To  obtain  an  explicit  formula  for  the  error,  lot  ns  assume  that  the  random  signal  has  a  fiat 
spectral  density: 


/>(") 


1  if  \v\  <  f 

0  if  M  >  f 


(17) 


With  a  flat  spectrum,  the  integrals  in  Equation  (l(i)  become  trivial  to  evaluate  analytically.  Com¬ 
bining  with  Equation  (13)  produces  the  following  asymptotic  expression  for  the  mean  square  error: 


(  ~ 


Nw 

N 


2p+\ 


Fn 


At 


P\A  t. 


ivq 


,  2/i At  1  as  N  — *  oo. 


where 


T 

N  =  — 

At 


Nw  = 


Et. 


2;> 4  I 


(2j>  f  l)7T2f>+2  I  w(t)'2(It 


FP(x,y)  =  -(Qp(x  +  y)  +  Qp(x-y)) 


Qp(z)  =  E 


1 


n  -oc 

7J  /  0 


,  (2«  -  z) 


-  xt2/>+1  ' 


(18) 


(19) 

(20) 
(21) 
(22) 


The  quantities  above  all  have  simple  interpretations.  N  is  the  number  of  samples  in  the 
window  and  is  not  required  to  be  an  integer.  The  choice  of  window  function  determines  the  order 
p  and  the  constant  Nu .  The  function  Fp  is  independent  of  the  window  function;  it  determines  the 
impact  of  the  normalized  sample  spacing  ^  and  scaled  frequency  //A/. 


5.1  EVALUATION  OF  QP(Z) 

For  p  >  0.  the  sum  in  Equation  (22)  is  absolutely  convergent  .  For  p  =  0.  the  terms  for  positive 
and  negative  n  may  be  combined,  yielding  a  convergent  series  which  is  evaluated  by  the  method  of 
residues: 


«»«=  £  5^TT  =  E 


n—  1 


2z  1  7T  HZ 

4n2  -  c2  _  7  ~  2  C  Ot  T' 


(23) 


Differentiating  the  above  identity  term-by-term  with  respect  to  z  leads  to  analytic  expressions 
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for  Qp.  p  >  0: 


«■(--)  =  ji  (t) 

<hU)  =  3  -  is(7r;j  +  5)  cot  (3)  l-srl  (  y  ) 

1  7r' (56 cos(ttz)  +OO.s(27Tc)  +  123) cot  (4f)  CSC*’  ) 


It  is  also  helpful  to  have  a  qualitative  understanding  of  the  behavior  of  Qp(z).  Due  to  the 
restrictions  At  <  A/My<(  and  |//|  <  •§ ,  it  is  clear  that  only  the  range  0  <  z  <  2  is  relevant.  From  the 
definition  in  Equation  (22),  it  is  easy  to  demonstrate  the  following  basic  properties  for  all  p  >  0: 


Qp(0)  =  » 

Qp(i)  =  1 

Jim  Qp{z)  =  oc 


Qp{s)  >  0  for  z  >  0 

Q'p{~)  >  ()  for  2  >  0 

(%(z)  =  (2p+l)(2p  +  2)Qr+l{z). 


It  directly  follows  that  all  derivatives  of  Qp(z)  are  positive  for  all  p  >  0.  2  >  0.  The  properties 
described  above  are  easily  seen  in  the  plots  of  Qp(z )  shown  in  Figure  1. 

For  2  <<  1,  the  analytic  forms  for  Qp{z)  given  above  become  difficult  to  compute  numerically 
due  to  the  cancellation  of  diverging  factors  as  c  — *  0.  Instead,  it  becomes  simpler  to  compute  using 
the  Taylor  series  around  2  =  0: 


Quiz) 


(p  +  hK(%P  +  2)  ,  jp  +  7}){p  +  FHp+  |)C i~P  +  4) 


4  p 


■2  + 


(j  -  4 p 


+  0(z% 


(24) 


where  C(.s)  = 


_  v^oc 


-1  n 


is  the  Rieniami  /eta  function. 


5.2  QUALITATIVE  DESCRIPTION  OF  ERROR 

The  analysis  above  loads  to  the  following  general  conclusions  regarding  the  asymptotic  error: 

1.  c2  decreases  like  1  ^  for  hxed  A/. 

2.  r2  decreases  nionotonically  with  At  for  fixed  N. 

3.  f2  increases  nionotonically  with  |//|. 

4.  f2  ~  1  for  fl  =  0  and  critically  sampled  data. 
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5.3  EVALUATION  OF  Nxy  FOR  STANDARD  WINDOW  FUNCTIONS 


Tabic  1  list  results  for  some  of  the  most,  common  window  functions.  For  each  window,  the 
order  p  and  constant  Nw  is  shown.  When  possible,  the  exact  analytic  expression  for  A>/H  3  is  also 
givcni.  In  some'  eases,  a  second  line  in  the  table  gives  the  next  order  term  for  the  same  window.  1  ho 
order  /;  is  determined  by  consecutively  examining  iy(r)  and  its  derivatives  until  a  discontinuity  is 
found.  The  constant  Nw  is  then  computed  directly  from  Equation  (20).  All  analytic  and  numeric 
results  were  generated  using  Mathematica  5.2. 

It  should  be  stressed  that  the  results  for  iiext-to-leading  order  terms  in  Table  1  are.  strict  ly 
speaking,  not  the  next  term  in  a  correct  asymptotic  expansion.  Instead,  they  are  simply  the  result 
obtained  by  pretending  that  the  leading  order  Nw  is  zero,  and  continuing  with  the  leading-order 
analysis  to  larger  p  until  a  non-zero  term  is  obtained.  This  is  useful  mainly  in  cases  where  the  true 
leading  order  Nw  is  numerically  very  small,  such  as  in  the  case'  of  the  Hamming  window. 


5.4  VERIFICATION  OF  ASYMPTOTIC  RESULTS 

Since  the  results  derived  an'  asymptotic  in  N.  there  is  no  guarantee  that  they  will  achieve 
any  specific  level  of  accuracy  except  in  the  limit  N  —  oc.  To  establish  the  utility  of  the  asymptotic* 
formula,  the  normalized  error  was  computed  via  Monte  Carlo  with  HP  random  trials.  Random  data 
was  generated  with  the  required  correlation  properties  at  a  sample'  spacing  of  ^y.  The  integral  in 
Equation  (6)  was  approximated  using  the  highly  oversampled  data.  The  sum  in  Equation  (7)  was 
computed  directly  using  every  20th  sample.  Averaging  over  values  of  o  was  approximated  by  a 
discrete  average  over  the  values  a  =  0.  . . .  ~.  The  Monte  Carlo  results  were  then  compared 

with  the  theoretical  asymptotic  predictions  from  Equation  (18). 

In  Figures  2  and  3?  the  errors  are  shown  as  a  function  of  number  of  samples  for  the  square, 
triangle.  Hanning.  Reisz,  and  Bolmian  windows.  For  these  examples,  the  frequency  is  //  =  0  and 
sampling  rate  is  Af/A/nvq  =  0.8.  The  agreement  between  Monte  Carlo  and  the  asymptotic  curves 
appears  to  be  excellent,  even  down  to  very  small  numbers  of  samples.  In  particular,  the  square 
window  results  match  Monte  Carlo  closely  even  for  a  single  sample. 

For  the  Hamming  and  Kaiser  windows,  the  leading  order  coefficient  Y„.  is  very  small.  Conse¬ 
quently,  convergence  to  the  asymptotic  result  requires  larger  values  of  A\  Figures  4  and  5  compare' 
Monte  Carlo  results  with  leading  order  asymptotics  as  well  as  leading  plus  subleading  order.  The 
subleadiug  term  appears  to  improve  the  results  for  small  A\  but  by  N  -  5  the  leading  term  suffices. 

In  Figure  (i,  errors  are  shown  as  a  function  of  frequency  for  the  square,  triangle,  and  Hanning 
windows.  The  number  of  samples  is  fixed  at  N  =  10  in  all  cases.  The  asymptotic  error  increases 
monotonically  with  |/i|.  as  described  earlier  in  Section  5.2.  While  the  agreement  is  good  near 
zero  frequency,  it  appears  that  the  asymptotic  formula  overestimates  the  error  near  the  Nyquist 
limit,  especially  for  the  triangle  and  Hanning  windows.  This  is  not  surprising,  since  the  asymptotic 
formula  ceases  to  be  valid  past  the  Nyquist  limit. 
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TABLE  1 


Parameters  for  Some  Standard  Window  Functions 


Name 

?e(r)  for  |t|  <  ^ 

P 

n'£p+1 

Nu, 

square 

i 

0 

2 

0.2026 

triangle 

1  -2|r| 

1 

24 

JF 

0.62G9 

Reisz 

1  -  At1 

1 

20 

0.5899 

Hamming 

0.54  +  0.4C  cos  (2ttt) 

0 

64 

10877T2 

0.0032G3 

2 

33856 

[mr,TT2 

0.8084 

Hanning 

0.5  4-  0.5  cos  (27rr ) 

2 

61 

157T2 

0.84 56 

Blackman 

0.42  +  0.5  cos  (2ttt)  + 
0.08  cos  (47tt) 

2 

5181 

76I5tt2 

0.5858 

Bohrnan 

(1  —  2|r|)  cos  (27 rr)  + 

^  sin  (2tt|t|) 

•? 

18432 

1.338 

♦  > 

77TJ(  15+27T2) 

Kaiser 

(4=0 

/o(v/1-'1t2) 

/o(l) 

0 

0.1460 

Kaiser 

(0=3) 

Kaiser 

(4=5) 

/0(W1-'It*) 

/o(3) 

/o(5V,1-4tj) 

/o(5) 

0 

1 

0 

1 

0.01599 

0.3523 

0.000677 

0.2426 
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Figure  I.  Plot  of  the  function 
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Figure  2.  Normalized  mean  squared  error'  versus  number  of  samples  in  window.  Results  are.  shown  for  both 
Monte  Carlo  and  the  asymptotic  error  form mla  in  Equation  (IS).  All  results  an  for  zero  frequency  (p  =  0 ) 
and  sampling  rate  of  At / At.nyq  =  0.8. 
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Fujun  l.  Normalized  mean  squared  envr  versus  number  of  samples  in  vrindow.  Results  an  shown  for  both 
MonU  Carlo  and,  the  asymptotic  error  foimnla  in  Equation  (IS).  All  ir. suits  an  for  zeiv  fnqueney  (fi  —  0) 
and  sampling  rate  of  At/ Atuyi]  —  0.8. 


N  (samples) 


Figure.  4 ■  Normalized  mean  squared  error  for  the  Hamming  window.  Results  are  shown  using  both  Monte 
Carlo  and  the  asymptotic  error  formula  in  Equation  (18).  Asymptotic,  results  air.  shown  for  both  leading 
order  and  leading  plus  subleading  order.  All  eases  acre,  for  zero  frequency  (ft  =  0)  and  sampling  rate  of 
At/A*nvq  =  0.8. 
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N  (samples) 


FiffUie  5.  Normalized  mean  squared  error  for  the  Kaiser  window  (R  5).  Results  tin  shown  using  both 
Monte  Carlo  and  the  asymptotic  error  formula  in  Equation  (IS).  Asymptotic  n  suits  an  shown  for  both 
leading  order  and  leading  plus  subleading  order.  All  eases  axe  for  zero  frequency  (p  0 )  and  sampling  rate 
ofAt'/Atu,„  -0.8. 
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Figuie  6.  Normalized  mean  squared  error  versus  frequency.  Results  are  shown  using  both  Monte  Carlo  and 
the  asymptotic  error  formula  in  Equation  (18).  AH  cases  are  for  N  =  10  samples  and  sampling  rate  of 
At./ At'nyq  =  0.8. 
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6.  CONCLUSIONS 


The  windowed  DFT  can  be  regarded  as  an  approximation  of  the  continuous  Fourier  transform 
with  the  same  window  function.  The  mean  square  error  for  stochastic  signals  is  for  the  most  part 
accurately  captured  by  the  asymptotic  formula  in  Equation  (18).  Parameters  for  some  of  the  most 
common  window  functions  are  tabulated  in  Table  1.  Monte  Carlo  calculations  confirm  that  the 
asymptotic  error  is  accurate,  even  when  the  number  of  samples  Ar  is  small.  Deviations  from  the 
asymptotic  result  mainly  occur  for  frequencies  approaching  the  Nyquist  limit. 

The*  asymptotic  error  formula  should  prove  especially  useful  in  scenarios  where  the'  processing 
parameters  must  he  determined  dynamically  during  processing.  In  such  case's,  the  asymptotic  result 
permits  one  to  analytically  optimize  any  design  metric. 
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APPENDIX  A 

MATHEMATICA  CODE  FOR  GENERATING  ASYMPTOTIC  RESULTS 


The  results  presented  in  Table  1  were  prepared  using  Mathematics  version  5.2.  Figure  A-l 
shows  the  actual  commands  used  to  generate  the  results,  along  with  comments  to  indicate  the 
purpose  of  various  steps.  Inputs  to  Mathematical  an'  indicated  in  bold  type,  while  outputs  and 
comments  an'  in  normal  type. 
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alias  1 3.  nb 


I 


Calculations  for  Asymptotic  Error  in  DFT 


Define  oversampling  error  function: 
q[p_*  *_]  :■ 

Sum[  1  /(2n-x)A(2p+l)-l/(2n  +  x)A(2pD),  (n,  1,  Infinity)]  //  FullSimplify 
Compute  first  few  cases: 
q[0,  x] 

i  i  ~  r  *  X  1 
x  2  nCO,-\  2  I 


q[i.  *1 


ST'  8  "  COtt  2 


Csc I 1 


q[2,  x] 


1  *  ,  r  ~  „  r  /T  X  i  „  r/tx 

X77  ’  T92  *  b  *  COS  1 CoL  [  ~2~  I  Csc  l  '  ■/ 


q[3  x] 

1  r’  { 123  ♦  56  Cos  |  n  x )  ♦  Cos  [  2  n  x]  )  Cot  ~  ]  Csc  ( 


x~  23040 

q[4,  x] 

1  <5786  +  404  7  Cos  [r,  xi  +  24  6  Cos|2  n  x]  ♦  Cos[3ti  x |  )  Cot  [  |  Csc  |  ^  |* 


5160960 


q[5,  x] 


x11  1857945600 

(tt11  (450995  *  408364  Cos|*  x]  *  46828  Cos(2  n  x]  ♦  1012Cos(3  n  x]  ♦  Cos[4  rr  xj  )  Cot.  [  ^  ' 


Define  window  functions  for  -1/2  <  x  <  1/2: 


sguara  [x_]  :  =  1 

hanning[x_]  :  =  1/2  f  (1/2)  Cos [2  •  Pi  *x] 

hamming [x_]  :=  (54  /  100)+  ( 46  /  100)  *  Cos  [2  *  Pi  *  x] 

blackman [x_]  ;=  42  /  100  +  (50  /  100)  Cos[2*Pi*x]  ♦  (8/100)  *  Cos  [4  *  Pi  *  x] 


Fifjuir.  A-l.  TmnscTvpi  of  Mathcviatica. 


22 


alias  li.nb 


aba  [x_]  :  =  Piecewise [  { (-x,  x<0},  {x,  xiO)}] 

.vn<  :  -i .  -f  t  1 .  V  »t  *p*»  .  i.n«i  «s»  r  r  -  rtvw  ,ir»  «  .  n.in-  "  »l  “  *  ;m.  \\  u*  »»..M  i  t  "At  “  M  rC... 

triangle [x_]  :  =  1  -  2  *  aba [x] 
reiaz[x_]  :  =  l-4xA2 

bohman  [x_]  :  *  (1-2  aba[x]  )  Coa  [2  Pi  x]  ♦  (1  /  Pi)  Sin  [2  Pi  aba  [x]  ] 
kaiaer[0_]  [x_]  :=  Beaaell[0,  0Sqrt[l-  (2  x)  A2]]  /Beaaell[0,  J9] 

Define  functions  to  compute  leading  order  error  term: 
diac[f_]  :=Module[(g,  x,  a}, 

g [y_]  * 

(Limi t[f  [x]  ,  x  -  >  y ,  Direction  ->  -  1]  -  Limit  [f  [x]  x  -  >  y  Direction  -  >  1  ]  )  /  /  Simplify 
a  *  Solve[Reduce [g [x]  !s  0,  x]  ,  x] ; 

Limit  [f  [x]  ,  x->l/2]A2  ♦ 

Limit  [f  [x]  x  ->  - 1  /  2]  A  2  ♦  If  (Length  [a]  >  0,  Total  [  (g[x]  /  .  a)  A2]  ,  0] 

1 

err[g0__]  :  *  Module  [  [x,  norm  g  ,  coef,  p  .  p2  .  temp  temp2) 
norms  Integrate [gO [x] A  2  ,  (x,  -1/2,  1/2}]; 
g[x_]  *  gO [x] ; 
coef  =  diac[g] ; 

p  =  0 ; 

While  [coef  si  0, 

g  [x_]  =  FullSimplify[D [g [x]  ,  x]  ,  x  €  Reala]  ;  coef  =  diac  [g ]  ;  p  =  p  ♦  1]  ; 
temp  =  Simplify!  (coef  /  (norm  *  (2p  +  1)  *  Pi  A  (2  *  p  ♦  2)  )  )  ]  ; 
g[x_]  *  FullSimplif y[D [g [x] ,  x]  x  €  Reala] ; 

{p,  temp,  N [  temp A  (1  /  (2  p  *  1)  )  ]  } 

] 

Define  function  to  compute  first  two  leading  error  terms 

err2[g0  ]  -Module[(x,  norm,  g,  coef,  p,  p2  ,  temp,  temp2)  , 
norm*  Integrate  [gO  [x]  A  2 ,  (x,  -1/2,  1/2}]; 
g[x_]  =  g0[x] ; 
coef  =  diac[g] 

P  *  0; 

While[coef  ss  0 , 

g  [x_]  *  FullSimplify[D  [g  [x]  ,  x]  ,  x  e  Reala]  ,  coef  «  diac  [g]  ;  p  *  p  +  1] 
temp  =  Simplify  [  (coef  /  (norm  *  (2  p  +  1}  *  Pi  A  (2  *  p  ♦  2) ) )  ] 
g[x_]  =  FullSiroplify [D [g [x]  x]  x  €  Reala] ; 
coef  =  diac[g]; 

p2  =  p  +  1 ; 

While  [coef  3=  0  , 

g[x_]  ■  FullSimplif  y[0  [g  [x]  ,  x]  ,  x  €  Reala]  coef  *  diac[g]  p2  =  p2  ♦  1] 
temp2  =  Simplify!  (coef  /  (norm  *  (2  p2  ♦  1)  *  Pi  A  (2  *  p2  ♦  2)  )  )  ]  ; 

{p,  temp,  N[tempA  (1  /  (2  p  ♦  1)  )  ]  ,  p2  ,  temp2  N[temp2A  (1  /  (2  p2  ♦  1)  )  ]  } 

] 


Figure  A- 1.  Transcript  of  Alatheinatica  (continued). 
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alias  13.  nh 


3 


Compute  results: 

err [square] 

JO,  0.202642} 


err2 [hanmng] 

\2,  rr^-#  0 . 845586,  4,  *-^4  ,  1.16134) 

1  1  5  /T*  2  7  n- 


err2 [hamming] 

Tsff^'  o-003263”'  -  °-80M1,i 


err2 [blackman] 


o  5184 
761b n' 


0.585786, 


173056 

15237^ 


i. 31193} 


err2 [reisz] 

jl,  ^  ,  0.589943,  2,  48  ,  0.549121} 

>  ;t 4  n*  1 

err [triangle] 

jl,  -r4  ,  0.626900) 


err2  [bohman] 

18432 

7  tt“  FT  5  ♦~T7^T ' 


1.33807,  5, 


1179648 
165  n-1  -.  22  H  ' 


1.68593} 


err2 [kaiser [1]  ] 


(0, 


_ _ _ 2 _ 

n* Bessel  I (0,  1  ] *’  J‘  —  tlx$  1 83  1  98 

_ 2_ _ 

3  7T4  Bessel  1 1  0,  IJ*T  ‘  cix$183198 


0. 146048, 


,  0.170226) 


err2 [kaiser [2] ] 

j0f  -  -  ,  0.0588426, 

7T-  Bes so  11(0,  2  ]  -  I  jr^  —  cl  x $  1  8 4  4  68 

1, - — — - - , - — - - ,  0.31680  9) 

3  Besse  1 1  [  0,  |  ’  * *'  *  — cl xSl  84468 


Figtif'c  A-l.  Transcript  of  Mat  hem  atica  (continued). 
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alius  1 3  nh 


err2 [kaiser [3] ] 


1 


0, 


Bes.se 3 1  0,  '| 


H  •  tlljl  . 


ft  4  »H*7  1 
j*ii  io".  f 


,  0 . 0  J  b  8  6  3 , 


•lx  .  8 


n<  1C,  i[- 


54 


. 


t"  •  ■  |  4  k;  IK  • 

*rpw»U|rt7;|' 


ci  xC>  1  8 1  \h 


err2 [kaiser [4 ] ] 

n  Be » sell 1 0.  4  —  '  ‘ ‘---4 --  .V**?7”*  ♦  dx$18 7 D © 4 


3  "4  Bessel  l  (  0 ,  4  | 3 


lilL^ 


4  J-  cl  x$  l  8^004 


err2 [kaiser [5] ] 

[0, 


n£  Bessel  I  [  0,  5  ' 


- ,  0 .0008  ' 

dx$  18  82 


3  n*  lk?S3C  1  ;  0  •  D|* 


*  »*»|o**  M  4»MlSJ7Jr  \S 

b-  -UK'  - 


3 .  ?  4  2  67  5  j 


tJx$l  B8?2." 


Figure  A-l.  Transcript  of  Mathematics  (continual). 
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